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The limit of neutron-rich nuclei, the neutron drip-line, evolves regularly from light to medium- 
mass nuclei except for a striking anomaly in the oxygen isotopes. This anomaly is not reproduced 
in shell-model calculations derived from microscopic two-nucleon forces. Here, we present the first 
microscopic explanation of the oxygen anomaly based on three-nucleon forces that have been estab- 
lished in few-body systems. This leads to repulsive contributions to the interactions among excess 
neutrons that change the location of the neutron drip-line from to the experimentally observed 
^*0. Since the mechanism is robust and general, our findings impact the prediction of the most 
neutron-rich nuclei and the synthesis of heavy elements in neutron-rich environments. 

PACS numbers: 21.10.-k, 21.30.-x, 21.60.Cs, 27.20.4-n 



One of the central challenges of nuclear physics is to 
develop a unified description of all nuclei created in the 
laboratory and the cosmos based on the underlying forces 
between neutrons and protons (nucleons). This involves 
understanding the sequences of isotopes in the nuclear 
chart. Fig. [TJ from the limits of proton-rich nuclei to 
the neutron drip-line. These limits have been established 
experimentally up to oxygen with proton number Z—S. 
Mapping out the neutron drip-line for larger Z ,1] and ex- 
ploring unexpected structures in neutron-rich nuclei are a 
current frontier in the physics of rare isotopes. The years 
of discovery in Fig. [T] highlight the tremendous advances 
made over the last decade. 

Figure [T] shows that the neutron drip-line evolves reg- 
ularly with increasing proton number, with an odd-even 
bound-unbound pattern due to neutron halos and pairing 
effects. The only known anomalous behavior is present 
in the oxygen isotopes, where the drip-line is strikingly 
close to the stability line ■ Already in the fluorine iso- 
topes, with one more proton, the drip-line is back to the 
regular trend Q. In this Letter, we discuss this puzzle 
and show that three-body forces are necessary to explain 
why ^''O [4, 5] is the heaviest oxygen isotope. 

Three-nucleon (3N) forces were introduced in the pio- 
neering work of Fujita and Miyazawa (FM) [3] and arise 
because nucleons are composite particles. The FM 3N 
mechanism is due to one nucleon virtually exciting a sec- 
ond nucleon to the A(1232MeV) resonance, which is de- 
excited by scattering off a third nucleon, see Fig. |3je). 

Three-nucleon interactions arise naturally in chiral ef- 
fective field theory (EFT) which provides a system- 
atic basis for nuclear forces, where nucleons interact via 
pion exchanges and shorter-range contact interactions. 
The resulting nuclear forces are organized in a system- 
atic expansion from leading to successively higher orders. 
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FIG. 1: Stable and unstable nuclei with Z ^ 14 and neutron 
number A'' 01 . The oxygen anomaly in the location of the 
neutron drip-line is highlighted. Element names and years of 
discovery of the most neutron-rich nuclei are given. The axis 
numbers indicate the conventional magic numbers. 



and include the A excitation as the dominant part of the 
leading 3N forces [8] . The quantitative role of 3N interac- 
tions has been highlighted in recent ab-initio calculations 
of light nuclei with A = N + Z ^^,12^^. 

We first discuss why the oxygen anomaly is not re- 
produced in shell-model calculations derived from micro- 
scopic NN forces. This can be understood starting from 
the stable '^^O and adding neutrons into single-particle 
orbitals (with standard quantum numbers nlj) above the 
^^O core. We will show that correlations do not change 
this intuitive picture. Starting from '^^O, neutrons first 
fill the orbitals, with a closed subshell configuration 
at {N = 14), then the lsi/2 orbitals at ^^O {N = 16), 
and finally the 0^3/2 orbitals at ^^O {N = 20). For sim- 
plicity, we will drop the n label in the following. 

In Fig. [21 we show the single-particle energies (SPE) of 



the neutron (1^/2, Si/2 and c?3/2 orbitals at subshell clo- 
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(c) G-matrix NN + 3N (A) forces 
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(d) V|^^^ ^ NN + 3N (A.n'lO) forces 
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FIG. 2: Single-particle energies of the neutron dr,/2, Si/2 and 
1^3/2 orbitals measured from the energy of '^^O as a function of 
neutron number A'^. (a) SPE calculated from a G matrix and 
from low-momentum interactions Viowfc- (b) SPE obtained 
from the phenomenological forces SDPF-M [T3| and USD- 
B [l^. (c,d) SPE including contributions from 3N forces due 
to A excitations and chiral EFT 3N interactions at N^LO [2^ . 
The changes due to 3N forces based on A excitations are 
highlighted by the shaded areas. 



sures N = 8, 14, 16, and 20. The evolution of the SPE 
is due to interactions as neutrons are added. For the 
SPE based on NN forces in Fig. [5] (a), the ^3/2 orbital 
decreases rapidly as neutrons occupy the 1^5/2 orbital, 
and remains well-bound from A'^ = 14 on. This leads 
to bound oxygen isotopes out to = 20 and puts the 
neutron drip- line incorrectly at ^*0. This result appears 
to depend only weakly on the renormalization method 
or the NN interaction used. We demonstrate thi s by 
showing SPE calculated in the G matrix formalism 



HI, 



which sums particle-particle ladders, and based on low- 
momentum interactions Mowfe 13 obtained from chiral 



NN interactions at next-to-next-to-next-to-leading order 
(N'^LO) [l^l using the renormalization group. Both cal- 
culations include core polarization effects perturbatively 
(including diagram Fig. |3] (d) with the A replaced by a 
nucleon and all other second-order diagrams) and start 
from empirical SPE lj| in ^^O. The empirical SPEs con- 
tain effects from the core and its excitations, including 
effects due to 3N forces. 

We next show in Fig. [5] (b) the SPE obtained from the 
phenomenological forces SDPF-M ^ and USD-B ^ 
that have been fit to reproduce experimental binding en- 



ergies and spectra. This shows a striking difference com- 
pared to Fig. [2] (a): As neutrons occupy the ^5/2 orbital, 
with A'^ evolving from 8 to 14, the ^3/2 orbital remains 
almost at the same energy and is not well-bound out to 
A^ = 20. The dominant differences between Figs. [2] (a) 
and (b) can be traced to the two-body monopole compo- 
nents, which determine the average interaction between 
two orbitals. The monopole components of a general two- 
body interaction V are given by an angular average over 
all possible orientations of the two nucleons in orbitals Ij 
and /'/ 0, 
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where the sum over magnetic quantum numbers m and 

[11 for 



can be restricted by antisymmetry (see [Ij 



details). The SPE of the orbital j is effectively shifted by 
y^rnono jnultiplied by the occupation number of the orbital 
' This leads to the change in the SPE and determines 
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shell structure and the location of the drip-line 

The comparison of Figs. [5] (a) and (b) suggests that the 
monopole interaction between the ^3/2 and ^5/2 orbitals 
obtained from NN theories is too attractive, and that the 
oxygen anomaly can be solved by additional repulsive 
contributions to the two-neutron monopole components, 
which approximately cancel the average NN attraction 
on the orbital. With extensive studies based on NN 
forces, it is unlikely that such a distinct property would 
have been missed, and it has been argued that 3N forces 



may be important for the monopole components [21 



Next, we show that 3N forces among two valence neu- 
trons and one nucleon in the ^^O core give rise to repul- 
sive monopole interactions between the valence neutrons. 
While the contributions of the FM 3N force to other 
quantities can be different, the shell-model configurations 
composed of valence neutrons probe the long-range parts 
of 3N forces. The repulsive nature of this 3N mechanism 
can be understood based on the Pauli exclusion princi- 
ple. Figure [3] (a) depicts the leading contribution to NN 
forces due to the excitation of a A, induced by the ex- 
change of pions with another nucleon. Because this is 
a second-order perturbation, its contribution to the en- 
ergy and to the two-neutron monopole components has 
to be attractive. This is part of the attractive ^3/2-^5/2 
monopole component obtained from NN forces. 

In nuclei, the process of Fig. |3] (a) leads to a change of 
the SPE of the j, m orbital due to the excitation of a core 
nucleon to a A, as illustrated in Fig. [3] (b) where the ini- 
tial valence neutron is virtually excited to another j', m' 
orbital. As discussed, this lowers the energy of the j, m 
orbital and thus increases its binding. However, in nuclei 
this process is forbidden by the Pauli exclusion princi- 
ple, if another neutron occupies the same orbital 
as shown in Fig. [3] (c) . The corresponding contribution 
must then be subtracted from the SPE change due to 
Fig. 13] (b). This is taken into account by the inclusion 
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FIG. 3: Processes involving 3N contributions. Tlie external 
lines are valence neutrons. The dashed and thick lines denote 
pions and A excitations, respectively. Nucleon-hole lines are 
indicated by downward arrows. The leading chiral 3N forces 
include the long-range two-pion-exchange parts, diagram (f), 
which take into account the excitation to a A and other res- 
onances, plus shorter-range one-pion exchange, diagram (g), 
and 3N contact interactions, diagram (h). 



of the exchange diagram. Fig. |3](d), where the neutrons 
in the intermediate state have been exchanged and this 
leads to the exchange of the final (or initial) orbital labels 
j,m and j',m'. Because this process reflects a cancella- 
tion of the lowering of the SPE, the contribution from 
Fig. [3] (d) has to be repulsive for two neutrons. Finally, 
we can rewrite Fig.[3](d) as the FM 3N force of Fig. [3] (e) , 
where the middle nucleon is summed over core nucleons. 
The importance of the cancellation between Figs. |3] (a) 
and (e) was recognized for nuclear matter in Ref. [22]. 

The process in Fig. [3] (d) corresponds to a two- valence- 
neutron monopole interaction, schematically illustrated 
in Fig. U (d). The resulting SPE evolution is shown in 
Fig. m (c) for the G matrix formalism, where a standard 
pion-N-A coupling [2^ was used and all 3N diagrams of 
the same order as Fig. [3] (d) are included. We observe 
that the repulsive FM 3N contributions become signifi- 
cant with increasing N and the resulting SPE structure 
is similar to that of phenomenological forces, where the 
orbital remains high. Next, we calculate the SPE 
from chiral low- momentum interactions Viowk, including 
the chaMes due to the leading (N^LO) 3N forces in chiral 



EFT 12J], see Figs. [3] (f)-(h). We consider also the SPE 
where 3N-force contributions are only due to A excita- 
tions [2^. The leading chiral 3N forces include the long- 
range two-pion-exchange part, Fig.[3](f), which takes into 
account the excitation to a A and other resonances, plus 
shorter-range 3N interactions. Figs. |3] (g) and (h), that 



range contributions due to A excitations dominate the 
changes in the SPE evolution and the effects of shorter- 
range 3N interactions are smaller. We point out that 
3N forces play a key role for the magic number iV = 14 
between ^5/2 and Si/2 27[, and that they enlarge the 
iV = 16 gap between si/2 and ^3/2 @. 

The contributions from Figs. [3] (f)-(h) (plus all ex- 
change terms) to the monopole components take into ac- 
count the normal-ordered two-body parts of 3N forces, 
where one of the nucleons is summed over all nucleons in 
the core. This is also motivated by recent coupled-cluster 
calculations [2^ , where residual 3N forces between three 
valence states were found to be small. In addition, the 
effects of 3N forces among three valence neutrons should 
be generally weaker due to the Pauli principle. 

Finally, we take into account many-body correlations 
by diagonalization in the valence space. The resulting 
ground-state energies of the oxygen isotopes are pre- 
sented in Fig. m Figure |4] (a) (based on phenomeno- 
logical forces) implies that many-body correlations do 
not change our picture developed from the SPE: The en- 
ergy decreases to = 16, but the ^3/2 neutrons added 
out to = 20 remain unbound. Figures S] (b) and (c) 
give the energies derived from NN forces, using a G ma- 
trix or low- momentum interactions Mow k , and including 
two- valence-neutron interactions due to 3N forces at the 
monopole level [2^. For all results based on NN forces, 
the energy decreases to = 20 and the neutron drip- 
line is incorrectly located at ^*0. The changes due to 3N 
forces based on A excitations are highlighted in Fig.|4](b) 
and (c). This leads to a better agreement with the ex- 
perimental energies and to a kink at iV = 16, which is 
further strengthened by shorter-range 3N forces, and for 
Fig.|4](c) leads to the neutron drip-line at ^"^O. 

The same 3N forces lead to repulsion in neutron mat- 
ter [s^]. Our results are also consistent with early shell- 
model explorations with 3N forces up to ^^O, where 
a small repulsive effect as in Figs. |4] (b) and (c) was 



have been constrained in few- nucleon systems 26|. The 



resulting SPE in Fig. [5] (d) demonstrate that the long- 



found [3]|. Because the formation of a halo is unre- 
alistic for the orbital and Si/2 is well bound (see 
Fig. [2] (b)), it seems unlikely that the ground states be- 
yond = 16 become bound by including the coupling 
to the continuum. This is consistent with Ref. [32|. We 
plan to study 3N-force effects on unbound states in the 
future using the methods of Refs. [H, Isl]. Fluorine iso- 
topes have one more proton than oxygen, and NN forces, 
primarily the tensor part, with this proton provide more 
binding to the valence neutrons [lo, Is^]. This valence 
proton-neutron effect is absent in the oxygen isotopes, 
making the repulsive 3N mechanism visible. Important 
directions for future work are to include the presented 3N 
contributions in coupled-cluster calculations [35| and in 
density-functional calculations, to systematically explore 
the effect over the full range of the nuclear chart. 

In summary, we have presented a robust 3N mecha- 
nism that provides repulsive monopole interactions be- 
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FIG. 4: Ground-state energies of oxygen isotopes measured from O, including experimental values of the bound O. 
Energies obtained from (a) phenomenological forces SDPF-M [TJ and USD-B [1^, (b) a G matrix and including FM 3N forces 
due to A excitations, and (c) from low-momentum interactions Viowk and including chiral EFT 3N interactions at N^LO as 
well as only due to A excitations [26l |. The changes due to 3N forces based on A excitations are highlighted by the shaded 
areas, (d) Schematic illustration of a two- valence-neutron interaction generated by 3N forces with a nucleon in the ^^O core. 



tween valence neutrons. Using microscopic NN and 3N 
forces as well as known SPE, our shell-model calculations 
naturally explain why ^^O is the heaviest oxygen isotope. 
The changes due to 3N forces are amplified and testable 
in neutron-rich nuclei and are expected to play a crucial 
role for matter at the extremes. 
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